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ABSTRACT 

The recent release of data on (i) high redshift source counts from NICMOS HUDF, and (ii) 
electron scattering optical depth from 3-year WMAP, require a re-examination of reionization 
scenarios. Using an improved self-consistent model, based on Choudhury & Ferrara (2005), 
we determine the range of reionization histories which can match a wide variety of data sets 
simultaneously. From this improved analysis we find that hydrogen reionization starts around 
z = 15, driven by the metal-free stars (with normal Salpeter-like IMF), and is 90% complete 
by z « 10. The photoionizing power of PopIII stars fades for z < 10 because of the con- 
comitant action of radiative and chemical feedbacks, which causes the reionization process 
to stretch considerably and to end only by z w 6. The combination of different data sets still 
favours a non-zero contribution from metal-free stars, forming with efficiencies > 2%. 
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1 INTRODUCTION 



The determination of the high Thomson electron scattering optical 
depth Tci = 0.17 ± 0.04 in the WMAP 1st year data (Kogut et al. 
2003; Spergel et al. 2003) had been a subject of extensive theo- 
retical study over the last few years. For sudden reionization mod- 
els, the high value of Td implied that reionization occurred at very 
high redshifts z ~ 15. This scenario seemed to be at tension with 
the measured Gunn-Peterson (GP) optical depth top (z ~ G) > 6 
from the absorption line experiments of 2 > 6 Sloan Digital Sky 
Survey (SDSS) QSOs (Fan et al. 2001; Fan et al. 2003; Fan et al. 
2005). Such high values of tgp seem to indicate that reionization 
was complete only at z ~ 6. 

Using a self-consistent formalism confronted with a wide 
range of observational data sets (redshift evolution of Lyman-limit 
absorption systems, GP and electron scattering optical depths, tem- 
perature of the mean intergalactic gas, and cosmic star formation 
history), Choudhury & Ferrara (2005; hereafter CF05) showed that 
the most favourable model is the one in which hydrogen reioniza- 
tion was complete at 2 ~ 12. By using the statistics of dark gaps 
in the Lya transmitted flux, this early reionization model was also 
shown not to be in conflict with QSO absorption line experiments at 
z > 6 (Gallerani, Choudhury, & Ferrara 2005). However, a nearly 
equally good fit to the data could be achieved for late reionization 
scenarios too (i.e., scenarios in which reionization is complete only 
at z ~ 6), particularly if one relaxes the constraints on Td. Given 
this, the need was to find an additional way to distinguish between 
the early and late reionization models, either through new theo- 
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retical methods (Gallerani, Choudhury, & Ferrara 2005) or from 
additional observational constraints. 

Fortunately, two new sets of data have been made available 
recently which could help in constraining the reionization history. 
First is the observations of high redshift sources in the NICMOS 
HUDF (Bouwens et al. 2005), where the analysis indicate that the 
number of sources at z ~ 10 should be three or fewer. This in- 
evitably rules out the occurrence of very massive (> SOOMq) 
stars (Salvaterra & Ferrara 2005). The second set of observations 
is the release of 3-year WMAP data (Hinshaw et al. 2006; Page 
et al. 2006; Spergel et al. 2006), which gives a lower value of 
Tci — 0.09 ± 0.03, thus questioning very early reionization sce- 
narios. 

Given these new data sets, it is important to find out the 
updated constraints on reionization using self-consistent models. 
More importantly, one has to address the issue as to which reion- 
ization histories are still viable and which sources are responsible 
for it. In this work, we extend the model of CF05 incorporating 
some additional physics (like chemical feedback) thus reducing the 
number of free parameters. We then confront the model with a wide 
variety of available data sets (including the two most recent ones 
above) with the aim of identifying a set of parameter values which 
can fit all the data sets simultaneously. 



2 MODEL DESCRIPTION 

In this Section, we first summarize the main features of the model 
introduced in CF05; following that we discuss the modifications 
and improvements made for the purpose of this work. 
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2.1 Summary of the model 

The main features of the semi-analytical model used in CF05 could 
be summarized along the following points (for detailed explana- 
tions see CF05). The model accounts for IGM inhomogeneities by 
adopting a lognormal distribution according to the method outlined 
in Miralda-Escude, Haehnelt, & Rees (2000); reionization is said 
to be complete once all the low-density regions (say, with overden- 
sities A < Acrit ~ 60) are ionized. Hence, the distribution of high 
density regions determines the mean free path of photons 

Ao 



Amfp(^) = 



(1) 



where Fy is the volume fraction of ionized regions and Ao is a nor- 
malization constant fixed by comparing with low redshift observa- 
tions. We follow the ionization and thermal histories of neutral, HII 
and Helll regions simultaneously and self-consistently, treating the 
IGM as a multi-phase medium. As reionization by UV sources is 
accompanied by heating of the meditmi, which can suppress star 
formation in low-mass halos, we compute the corresponding radia- 
tive feedback self-consistently from the evolution of the thermal 
properties of the IGM. To calculate the ionizing flux, three sources 
have been assumed: (i) metal-free (i.e., PopIII) stars, which dom- 
inate the flux at high redshifts. In CF05, they were asstuned to be 
massive (M > IOOMq); (ii) PopII stars with sub-solar metallic- 
ities having a Salpeter IMF in the mass range 1 — lOOM©. The 
transition redshift, «trans from PopIII to PopII stars was a free pa- 
rameter, usually fixed to be «trans ^ 10; (iii) QSOs, which are 
significant source for hard photons at 2 < 6; however they have no 
effect on the IGM at higher redshifts. 



2.2 New features 

We now discuss the additional physics we have incorporated in this 
work in order to further improve the predictive power of the model. 

• Radiative feedback: We have assumed that no haloes with 
virial temperatures lower than lO"* K are able to form stars; this 
completely neglects the contribution of minihaloes, which is now 
strongly supported by the 3 -year WMAP data (Haiman & Bryan 
2006). 

• Chemical feedback: The main limitation of CF05 model was 
the ideaUzed PopIII — » PopII transition which was assumed to start 
at ztrans and last for a dynamical time of the halo. According to the 
standard chemical feedback interpretation (Schneider et al. 2002; 
Schneider et al. 2003), the transition is driven by the enrichment 
of the medium which forces a drastic change in the fragmentation 
properties of star-forming clouds when metallicity exceeds the crit- 
ical value of Zcrit = 10^''^^ .^o (Schneider et al. 2002; Schneider 
et al. 2003). Such feedback-regulated transition has been studied in 
detail by Schneider et al. (2005), using a merger-tree approach to 
determine the termination of PopIII star formation in a given halo. 
We incorporate the same prescription in our model (using Fig 3 of 
Schneider et al. 2005), which allows us to compute the transition in 
a self-consistent manner. The main difference with respect to CF05 
is that the transition occurs over a prolonged epoch, i.e., no precise 
transition redshift can be identified. 

• IMF of PopIII stars: In CF05, a top-heavy IMF for PopIII stars 
was used, which was found to be disfavoured by a combination of 
constraints from source counts at 2; « 10 and the first year WMAP 
data (Schneider et al. 2005). In this work we use a very "conser- 
vative" assumption that the metal-free PopIII stars have a simple 



Salpeter IMF, just like the PopII stars, which is similar to the hy- 
pothesis made in Ciardi, Ferrara, & White (2003). One should keep 
in irund that the recent 3-year WMAP data need not necessarily 
rule out the possibility that PopIII stars have a top-heavy IMF; how- 
ever, wc limit ourselves to the most conservative model and check 
whether it can match all available observations. 

• Escape fraction: In CF05, the escape fractions for PopII and 
Popni stars, /esc.ii and /esc,iii, were considered as free parameters, 
independent of the halo mass, M, and redshift z. In reality, the 
situation is quite complex and the escape fractions do depend on 
both M and z. Unfortunately, there is still no good understanding 
of the process so as to model it theoretically. 

In this work, we retain the assumption that the escape fraction is 
independent of M and z; however, we use a physical argument to 
relate the escape fraction of PopII and PopIII stars. This is based 
on the fact that the escape fraction should scale according to the 
number of ionizing photons produced by a given source. Let A'abs 
denote that number of photons that can be potentially absorbed by 
the star-forming halo (which can be quite different from the num- 
ber of photons actually absorbed). It is usually proportional to the 
quantity CaR{T)nHne, where C is the clumping factor of the halo 
gas density inhomogeneities. There are further uncertainties related 
to the distribution of stars within the halo, and wc assume that such 
uncertainties can be absorbed within the proportionality factor. Let 
^■y,u (N^jii) denote the number of photons produced by PopII 
(PopIII) stars per unit mass of star formed and eN.,ii(e*,iii) denote 
the star-forming efficiency of the population. We can then define 
the parameter 



Ve: 



e*,iiAf-Y II 



(2) 



which measures the fraction of photons absorbed in the halo. Then 
one can write the relation 



fe: 



Min 



A^abs 



e*,iiN. 



7,11 



Min[l,77ei 



(3) 



which takes into account the fact that /esc, 11 = if j^osc > 1 (which 
essentially means that the halo is capable of absorbing more pho- 
tons than what is produced by the stars and thus all the photons 
produced are absorbed within the halo). Note that a higher value of 
e*,ii would give a higher /esc, 11 signifying that a higher fraction of 
photons will escape if the number of photons produced is larger. 

We now make the simplifying assumption that A'^abs depends 
only on the properties of the halo and is independent of the na- 
ture of the stellar source. Then the escape fraction for Popin stars 
would be given by 



1 - Min 

1 - Min 



iVabs 



e*,iiiiVT.,iii 



e*,iiiVT,,ii 

^ 1 ~ ^/csc 



e*,iiiA^. 



7,111 



(4) 

(5) 



which relates the escape fractions of the two stellar populations. 
Note that no assumptions about the gas density structure has been 
made; we simply used the fact that a higher fraction of photons will 
escape if the number of photons produced is larger. The above pre- 
scription can be extended to helium too. It thus helps us in reducing 
the number of free parameters in our model with the escape fraction 
being given by a single free parameter rjesc- 

• Self-consistent calculation of the temperature-density rela- 
tion: For calculations of the transmitted flux of the IGM (as would 
be observed in QSO absorption line experiments), it is usually as- 
simied that the temperature-density relation follows a power-law 
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form, T oc A''"^. In this work, we solve the temperature evolution 
equation for fluid elements of different densities and thus obtain 
the value of 7 at each redshift in a self-consistent manner (Hui & 
Gnedin 1997). 

• Additional observational constraints: In addition to the obser- 
vations described in CF05, we use a few additional constraints to 
determine our free parameters. The most notable of these is the ex- 
periments related to the source counts at high redshifts (Bouwens 
et al. 2005). Three possible high redshift candidates have been iden- 
tified by applying the I-dropout technique to the NICMOS HUDF; 
however the precise nature of these three sources could not be con- 
firmed. Hence, Bouwens et al. (2005) concluded that the actual 
number of z « 10 sources in the NICMOS parallel fields must 
be three or fewer. 

The number of sources above a redshift z observed within a solid 
angle dSl in the flux range [Fi^q , Fu^ + dF^^] is 

^-''o(>^) = dndK;(^-'^) = / ^^di^d^U;^^-'^)^^) 

where <W jdz dQ. denotes the comoving volume element per unit 
redshift per unit solid angle, and 

m:}^"-'^=\, dJi:(^-'^^'-*-)dMd^(^'^ 

is the comoving number of objects at redshift z' with observed flux 
within [F^Q , F^^ ^dF^^. The quantity d^njdMdz" gives the for- 
mation rate of haloes of mass M, calculated using Press-Schechter 
formalism. The flux is related to the mass of the halo M by the 
relation 

P _ e.(n,/»^)M/di/' -t./Q e-^eff(--o.^=o.^O 

47rdi(0')A!/o 

where e, is the star-forming efficiency of the population under con- 
sideration, 1^1 [tz' — t^ii ) is template luminosity per unit solar mass 
for the stellar population of age t^i — t^n (the time elapsed between 
the two redshifts), dL(z') is the luminosity distance and Az^o is the 
instrumental bandwidth. The quantity Teff(i-'o, z = Q, z') is the ef- 
fective optical depth at i^o between z' and z = 0, which can be 
calculated self-consistently from the semi-analytical model given 
the density distribution. While calculating the source distribution, 
we apply the same selection criteria as is used in the observational 
analysis. For calculating the template luminosity L, we use stellar 
population models of Bruzual & Chariot (2003) for Popn stars and 
of Schaerer (2002) for PopIU stars. 

We have also incorporated constraints from the observed trans- 
mitted flux in the Ly/3 region of the (}SO absorption spectra (in ad- 
dition to Lya), setting more severe constraints on the background 
ionizing flux. 

2.3 Free parameters 

For the background cosmology, we use the best-fit parameters as 
given by the 3-year WMAP data (Spergel et al. 2006), i.e., we as- 
sume a flat universe with total matter, vacuum, and baryonic densi- 
ties in units of the critical density of Q,m = 0.24, JIa = 0.76, and 
Qbh^ = 0.022, respectively, and a Hubble constant of Ho = 100 h 
km s~^ Mpc~^, with h = 0.73. The parameters defining the 
linear dark matter power spectrum are as = 0.74, Us = 0.95, 
dns/dlnk — 0. The reionization model, after improvements and 
modifications, contain just four free parameters. They are the star- 
forming efficiencies of PopII and Pom stars (e*,ii and e»,iii respec- 
tively), Tjesc related to the escape fractions of the two stellar pop- 



ulations [see equations (2)-(5)] and the normalization of the mean 
free path Ao [equation (1)]. The main exercise of this paper would 
be to find the range of parameter values which can match all the 
observational data we have considered. A detailed exploration of 
the parameter space using a statistical approach would be reported 
elsewhere. 



3 RESULTS 

We now constrain the above free parameters and select the best-fit 
model as the one that fits simultaneously all the available experi- 
mental data. 

3.1 The best-fit model 

The best-fit model is characterized by the parameter values e*,:: = 
0.2,e,,iii = 0.07,/csc,ii = 0.003. /e.cJii = 0.72 (keeping in 
mind that /esc.ii and /esc,iii are not independent). The comparison 
between the best-fit model and different observations is shown in 
nine panels of Figure 1. Quite remarkably, the model matches a 
wide variety of observations by fitting only four parameters. 

It is clear from the evolution of the neutral hydrogen frac- 
tion xm [Panel (a)] that current observations favour a model where 
reionization starts aroimd z = 15 and is 90 per cent complete by 
z ^ 10. 

The initial phase of reionization is driven by metal-free 
(PopIII) stars which are capable of producing a large number of 
ionizing photons, as can be seen from the evolution of the pho- 
toionization rate of neutral hydrogen Fp" [Panel (h)]. The Popin 
star-formation is severely quenched at z < 10 both because of the 
concomitant action of both radiative and chemical feedback. In fact, 
chemical feedback allows PopIII star formation only in small mass 
(< 2 X 10® M0) haloes, while radiative feedback tends to prohibit 
star formation in those, thus gradually terminating the PopIII stars. 
The total mass of PopIII stars formed is 0,jii ~ 10^®. As PopIII 
stars can no longer be formed, the progress of ionization fronts is 
limited, and hence the reionization extends and is completed only 
by 2 « 6. This evolution of xm is consistent both with high GP 
optical depths for Lya [Panel (e)] and Ly/3 [Panel (f)] and with the 
constraints from Lya emitters at z ~ 6.5 [Panel (a)]. 

At lower redshifts, the ionizing background seems to be dom- 
inated by QSOs as the Popn stars have negligible escape fraction. 
Unlike CF05, the escape fractions of the PopII and Popin stars in 
this work are determined by the same quantity ?7esc ; hence it is not 
possible to hike up /esc, 11 without affecting /esc.iii, in turn possibly 
violating some other observational constraints. At « w 3 — 5, our 
best-fit model is consistent with the observed evolution of the SFR 
[Panel (b)], Lyman-limit systems [Panel (g)] and the temperature 
of the mean density IGM [Panel (i)] ; finally, our estimates of Fpp 
are consistent with those obtained from numerical simulations of 
(Bolton et al. 2005) [Panel (h)]. Note that the SFR is almost always 
dominated by PopII stars [Panel (b)] as metal-free stars have a low 
star formation efficiency (3 times smaller than that of PopII stars); 
however, PopIII stars can still dominate the photoionization rate at 
high redshifts as they produce larger niunber of photons per imit 
mass of stars formed. 

The best-fit model produces negligible source counts at z » 
10, which suggests that the three NICMOS HUDF candidates at 
z « 10 are probably spurious detections. Note that the main con- 
tribution to the source counts comes from the nebular and Lya 
line emission of the Popni stars. The low value of source cotmt 
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Figure 1. Comparison of model predictions with obsei'vations for the best-fit model with parameter values e*,ii = 0.2, e, ni = 0.07, /esc, II = 
0.003, /osc, III = 0.72 (keeping in mind that /esc, II and /esc, III are not independent). The different panels indicate: (a) The volume-averaged neutral 
hydrogen fraction xpji, with observational lower limit from QSO absorption lines aX z = 6 and upper limit from Lya emitters at z = 6.5 (shown with arrows). 
In addition, the ionized fraction Xe is shown by the dashed line, (b) SFR for different stellar populations. The points with error-bars indicate low-redshift 
observations taken from the compilation of Nagamine et al. (2004) (c) The number of source counts above a given redshift, with the observational upper 
limit from NICMOS HUDF shown by the aiTow. The contiibution to the source count is zero at low redshifts because of the J-dropout selection criterion, (d) 
Electron scattering optical depth, with observational constraint from WMAP 3-year data release, (e) Lyo effective optical depth, with observed data points 
from Songaila (2004). (f) Ly/3 effective optical depth, with observed data points from Songaila (2004). (g) Evolution of Lyman-limit systems, with observed 
data points from Stonie-Lombardi et al. (1994). (h) Photoionization rates for hydrogen, with estimates from numerical simulations (shown by points with 
error-bars; Bolton et al. 2005). (i) Temperature of the mean density IGM, with observational estimates from Schaye et al. (1999). 



is mainly determined by the required high value of /csc,iii; since 
most of the photons escape the host halo, the amount of nebular 
and Lya line emission is small, and hence no sources are above 
the detection threshold of the NICMOS experiments. In this sense, 
one can rule out this best-fit model if at least one of the sources 
observed in the NICMOS HUDF turns out be indeed at z ^ 10. 
We shall discuss this possibility later. Before addressing other is- 
sues, it might be worthwhile mentioning the ionization history of 
doubly-ionized helium. We find that the escape fraction for photons 
with energies above 54.4 eV is not very high for PopIII stars, and 
hence the propagation of doubly-ionized helium fronts is not very 
efficient at high redshifts. The complete reionization occurs only 
around z ~ 3.5 because of QSOs. 



3.2 Variants of the best-fit model 

In spite of the success of our best-fit model in fitting observations, 
it is instructive to study some of its variants. In particular, we ask 
some interesting questions and try to find what the current data 
imply: 



(i) Is it possible to fit the data with reionization at higher red- 
shifts? By increasing e,,iii one can force an earlier start of the 
reionization process. For example, a model with parameter values 
e*,ii = 0.2, e,,ni = 0.2, /csc,ii = 0.006, /esc,iii = 0.9 gives a 
nearly equal good fit to the data. In this model, reionization starts 
much earlier and the IGM is 95 (99) per cent ionized by z ~ 10(8). 
However, as in the best-fit model, the contribution of the PopIII 
stars to the ionizing flux decreases because of feedback and hence 
the reionization is extended till 2 ~ 6. Interestingly, the most se- 
vere constraint on the early reionization scenarios comes from Ly/3 
observations at z ~ 6 which rule out high values of e*,iii,/csc,iii. 

(ii) What if one or some of the candidates in the NICMOS 
HUDF do turn out to be valid z ~ 10 sources? Then the best- 
fit model above could be ruled out as it predicts negligible source 
counts at z > 10. However, there are parameters within 2a of the 
best-fit value which predict high number of sources at 2 ~ 10. 
For example, a model with parameter values e»,n — 0.2, e,.iii — 
0.4, /esc, II = 0.0, /esc, III = 0.36 predicts « 2.6 sources at 
z ~ 10. The condition for producing high source counts is that 
the escape fraction of PopIII stars should be low, which in turn 
means that the ionizing flux is lower and reionization would be de- 
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layed. Such models tend to overpredict the Lya optical depth at 
« « 4 — 5 and thus the parameter space is severely constrained. In 
other words, if one observes sources at high redshifts, a very small 
parameter space would be allowed, and we would possibly be able 
to identify uniquely how reionization occurred. 

(iii) Do we still require a prolonged epoch of metal-free star 
formation to explain observations? In fact, one can explain the low 
redshift SFR constraints and the WMAP Tei without Popin stars 
provided e,,ii = 0.2 and /esc > 0.07. However, such a high 
value of escape fraction violates the Lya and Ly/3 optical depths 
at « K, 3 — 4; the GP optical depth measurements require that 
/esc < 0.05, which when combined with the WMAP t^i, gives 

> 0.02. Thus a combination of low redshift SFR, t^i and GP 
optical depth constraints imply that we still do require a non-zero 
contribution from metal-free PopIII stars, albeit with a small star- 
forming efficiency. 



4 SUMMARY 

We have extended the self-consistent model of CF05 incorporating 
key additional physical processes and compared the model predic- 
tions with a variety of data sets. Our formalism now includes the 
inhomogeneous IGM density distribution, three different classes of 
ionizing photon sources (metal-free PopIII stars, PopII stars and 
QSOs), chemical feedback inhibiting formation of PopIII stars in 
metal-enriched haloes, and radiative feedback preventing the for- 
mation of stars in galaxies below a certain circular velocity thresh- 
old. Our model is able to: (i) predict the star formation/emissivity 
history of sources and the number of sources above a given flux 
detection threshold at various redshifts, (ii) follow the evolution of 
H and He reionization and of the intergalactic gas temperature, and 
(iii) yield a number of additional predictions involving directly ob- 
servable quantities. 

By constraining the model free parameters with the available 
experimental data we have found a best-fit model and a set of al- 
lowed parameter values which matches very well all the available 
observations. From this analysis, an updated reionization scenario, 
which also takes into account the 3-year WMAP data, emerges: 

• Hydrogen reionization starts around z « 15 driven by metal- 
free (PopIII) stars, and it is 90 per cent complete by 2 ~ 10. 
The contribution of PopIII stars decrease below this redshift be- 
cause of the combined action of radiative and chemical feedback. 
As a result, reionization is extended considerably completing only 
at 2 ~ 6. 

• Scenarios in which reionization is completed much earlier are 
ruled out by a combination of constraints from Tei, the NICMOS 
source counts at 2 « 10 and the Ly/3 optical depths at 2 « 6. 

• The combination of Tei and GP optical depth constraints re- 
quire non-zero contribution from metal-free stars with a normal 
Salpeter IMF. Non-inclusion of PopIII stars would require a rel- 
atively larger ionization flux from PopII stars to match the WMAP 
Tel, which would then violate the GP optical depth constraints. 
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